Weber blockade theory of magnetoresistance oscillations in superconducting strips 
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Recent experiments on the conductance of thin, narrow superconducting strips have found periodic 
fluctuations, as a function of the perpendicular magnetic field, with a period corresponding to 
approximately two flux quanta per strip area [A. Johansson et al., Phys. Rev. Lett. 95, lf6805 
(2005)]. We argue that the low-energy degrees of freedom responsible for dissipation correspond 
to vortex motion. Using vortex/charge duality, we show that the superconducting strip behaves as 
the dual of a quantum dot, with the vortices, magnetic field, and bias current respectively playing 
the roles of the electrons, gate voltage and source-drain voltage. In the bias-current vs. magnetic- 
field plane, the strip conductance displays what we term 'Weber blockade' diamonds, with vortex 
conductance maxima (i.e., electrical resistance maxima) that, at small bias-currents, correspond to 
the fields at which strip states of N and TV + 1 vortices have equal energy. 



Introduction - It is often effective to characterize 
strongly correlated quantum systems in terms of the 
emergent, collective freedoms that describe their low- 
energy behavior. Vortices in superconductors consti- 
tute the most prominent example of such freedoms, and 
it has proven useful to address the Kosterlitz-Thouless 
phase transition exhibited by thin-film superconductors 
in terms of the statistical mechanics of an interacting 
plasma of such vortices. The supcrfluid-insulator quan- 
tum phase transition — exhibited by quite a number of 
systems [IMS), particularly granular InO films [U [S]— 
is widely suspected to result from a vortex prolifera- 
tion transition. If so, this behavior would provide a 
convincing demonstration of vortices behaving quantum- 
mechanically: not only do they exhibit quantal motion 
as individuals [5] , but they are also able to Bose- Einstein 
condense. 

Motivated by recent experiments by the Shahar group 
on superconducting InO strips [7], which exhibited os- 
cillations in the resistance as a function of magnetic 
field, in this paper we address such oscillations from 
the vortex point of view, and show that the notion of a 
"vortex blockade" allows us to explain such oscillations, 
in analogy with the Coulomb blockade theory that has 
been applied widely to electron passage through quan- 
tum dots. Our approach will be applied to the case 
of a vortex blockade in strips that are narrow. By us- 
ing charge-vortex duality, together with energetics argu- 
ments, we show that a vortex blockade results in a series 
of 'Weber' diamonds, which are analogs of the Coulomb- 
blockade diamonds, except for the important distinction 
that the electrical conductivity is maximal (rather than 
minimal) inside the Weber diamonds. We then apply 
the Beenakker formalism, developed to account for dis- 
sipation in quantum-dot transport [5], to the case of a 
superconducting strip, with the dissipation taking place 
via the normal modes of the vortex 'crystal' in the strip. 
Within this dissipative model, we predict that the resis- 
tivity of the strip should diverge — possibly observably — 



at low temperatures, as a power law. We conclude with a 
comparison between our results and experimental data, 
and thus show that our energetics considerations cor- 
rectly account for period of the magnetoresistance os- 
cillations, observed experimentally. 
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FIG. 1: (a) Schematic representation of the experimental 
setup. A superconducting strip (gray) is immersed in a mag- 
netic field B pointing out of the plane of the strip. The strip 
is contacted via a pair of leads (blue) that are used to pass a 
bias supercurrent J along the length of the strip. Simultane- 
ously, the time-averaged voltage is measured using the same 
leads. Vortices (depicted by swirls) reside in the strip, but 
occasionally traverse the strip (cf. the vortex trajectory indi- 
cated by arrows). Via the Josephson relation, such crossings 
correspond to voltage spikes, and thus result in dissipation, 
(b) Analogous quantum-dot circuit. The dot is depicted by 
the gray disc. The vortices and the magnetic field in the strip 
correspond, respectively, to the electrons and the gate voltage 
V g of the dot, the latter controlling the mean number of parti- 
cles on the dot. The source and drain for the vortices are the 
vacua adjacent to the strip. These vacua can be regarded as 
vortex condensates (green), with the potential energy differ- 
ence between the condensates set by the bias current. Thus, J 
in the vortex analogy corresponds to the source-drain voltage 
Vs — Vd in the dot (see text for details) . 
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Setup - We are concerned with describing the current- 
voltage characteristics of a superconducting strip that is 
subject to a perpendicular magnetic field. The setup that 
we have in mind is depicted in Fig. [IJi. It consists of a su- 
perconducting strip, contacted via a pair of leads that are 
used to inject a supercurrent into the strip. The voltage 
across the strip is read out using the same leads, and we 
rely on the assumption that the dissipation being probed 
originates entirely within the superconducting strip and 
not in the leads. 

We make several restrictions on the strip geometry that 
are pertinent from the vortex perspective. First, we treat 
the superconductivity in the strip as being essentially two 
dimensional, i.e., the strip thickness t is smaller than the 
coherence length £. This restriction implies that vortices 
in the strip are pancake (i.e., Pearl) vortices, as opposed 
to line vortices [5]. Vortices in thin strips have been in- 
vestigated previously, using magnetic force imaging |10j . 
These investigations reveal that strips are of two types: 
narrow and wide. In narrow strips (i.e., those having 
a width no greater than several coherence lengths), the 
vortices always lie in a line along the strip center, like 
peas in a pod. For increased magnetic fields, the linear 
density of the vortices is increased, until an upper crit- 
ical field is reached. On the other hand, in wider strips 
the vortices tend to form Abrikosov-like lattices at higher 
fields. In the present work, we concentrate on the case 
of narrow strips, and use intuition gained from Ref. [10] 
to describe the vortices. We shall address the question 
of wider strips in a separate article. Finally, we com- 
ment that we consider the strip to be free of sites that 
would pin vortices, and also to be of finite length (in 
practice [7J, ~ 10 — 100 £) so that configurations having 
distinct numbers of vortices have sufficiently different en- 
ergies as compared to fc^T. 

Dual description of a superconducting strip - Via 
charge/vortex duality, we map the behavior of a super- 
conducting strip to that of a quantum dot capable of 
undergoing a Coulomb blockade (see Fig. ]J>). In this 
mapping, the superconducting strip plays the role of the 
quantum dot. The vacua running along either side of 
the strip act as reservoirs that can exchange vortices 
with the strip, and thus play the role of the source and 
drain contacts of a quantum dot. As the vacua do not a 
have well-defined number of vortices, we describe them 
as vortex condensates of well-defined phase. The mag- 
netic field controls the chemical potential of the vortices 
on the strip, and thus acts as the gate voltage of the dot. 
Lastly, the bias current that runs along the strip exerts a 
Magnus force on the vortices, which acts in the direction 
perpendicular to the current flow. Effectively, this force 
creates a potential difference between the vacua on ei- 
ther side of the strip, and thus it acts as the source-drain 
bias of the dot does. The transport of vortices across the 
superconducting strip has a direct consequence for elec- 
trical transport within the strip, viz., each time a vortex 
passes across the strip, the phase between the left and 
right leads changes by 2ir. The voltage V between the 



leads is given by the Josephson relation: 
V = <P IV, 



(1) 



where Tn is the net rate of vortex tunneling, and $o = 
h/2e is the superconducting flux quantum. 

Having defined the analogy, we split the problem of 
computing the vortex conductance (i.e., the electrical re- 
sistance) of the strip into two sub-problems. The first 
consists of describing the energetics of static vortices. 
That is to say, we describe the energy cost of adding an 
extra vortex to a line of vortices already in the strip. 
The second sub-problem consists of utilizing the vor- 
tex energetics to construct a master equation describ- 
ing vortex hopping. In constructing the master equa- 
tion we make the typical — and experimentally relevant — 
assumption that vortex tunneling rates T t (Tb) between 
the strip and the top (bottom) vortex reservoirs are the 
smallest energy scales in the system, i.e., r t ,r;, <C ksT. 
This assumption ensures that after a vortex tunnels onto 
the strip it is completely de-phased, before it is likely to 
tunnel again, therefore validating the use of the master 
equation approach. 

Vortex energetics - To use the aforementioned analogy, 
we need to quantify the vortex energetics. To do this, we 
describe the strip in terms of a phase-only model Hamil- 
tonian, 



H = 



P 







j dx dy 





(2) 



in which p ~ fcs^BKT is the superfluid stiffness, <f> is the 
phase of the order parameter, and A = —B(y — w/2)e x 
is the vector potential in the London gauge. We use a 
coordinate system in which y runs from to w across 
the width of the strip and the long direction of the strip 
runs in the x direction, as indicated in Fig. [l^,. We sup- 
plement the Hamiltonian with the boundary conditions 
of no supercurrent through the top or bottom edges of 
the strip. The boundary conditions for the left and right 
ends of the strip depend on the properties of the contacts, 
and we return to this question, later. For the purpose 
of quantitative energetics, we invoke periodic boundary 
conditions in the long direction. We note, however, that 
any choice of boundary conditions that preserves the dis- 
creteness of the vortex "charging" energy would result in 
magnetoresistance oscillations. To proceed, we break the 
first sub-problem up into a single-vortex energetics part 
and an inter- vortex interaction part. 

The potential energy of a single vortex in a strip, as a 
function of its position y across the strip, was thoroughly 
explored previously (TT], with a result that comprises 
four parts, as follows: (i) The vortex core energy E cole , 
which is determined via the microscopic theory of the 
vortex, (ii) The energy of interaction of the vortex with 
its images, i.e., E- lma ,g e = — 27rpm[sin(7r£/ , u;)/ sm(iry/w)]. 
This interaction diverges logarithmically, as the vortex 
approaches the strip edge, due to the closest image- 
vortex; however, the divergence is unphysical, as vor- 
tices are finite in size, and one should therefore cut the 
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FIG. 2: (a) Vortex potential energy as a function of position 
across the strip for various magnetic fields for a 94 nm wide 
strip with £ = 18 nm and p — 5K, (see Ref. (b) Same 

as (a), but with nonzero bias current applied along the strip 
and magnetic field fixed at 0.12 T. 



divergence off at the coherence lengthscale £. In fur- 
ther calculations, we shall absorb the vortex core energy 
into the cut-off lengthscale. (iii) The part that accounts 
for the interaction of the vortex with the magnetic field, 



i.e. 



E 



vB 



= (2^) 2 (pB/$ ) (y-w/2Y-{Z-w/2y 



(iv) Finally, the part that describes the potential due 
to the Magnus force, i.e., -^Magnus = $o(J/w)(y - w/2), 
where J is the bias current. 

The consequences of these four parts are summarized 
in Fig. [2^,, in which we plot the potential energy of a 
single vortex in an infinitely long strip, as a function of 
its position. For small applied magnetic fields, vortices 
are unstable everywhere inside the strip. For larger mag- 
netic fields, the configuration having a vortex trapped in 
the middle of the strip becomes metastable. At a still 
larger field (i.e., the lower critical field) the vortex be- 
comes stable, globally. A bias current along the strip 
produces a tilting of the vortex potential energy diagram 
(see Fig. [2}}). In particular, we see that a vortex at the 
top edge of the strip has a potential energy that differs 
from that of a vortex at the bottom edge. Wc identify this 
energy difference, which corresponds to the work done by 
the Magnus force as a vortex is moved across the strip, as 
the difference in potential between the top and bottom 
vortex reservoirs. 

To find the interaction energy between vortices, de- 
scribed by Eq. (pi), we use the conformal mapping z — x+ 
iy = cxp(iru/w)to identify the strip < Imti < w and 
the entire upper half-plane < Imz < oo. The mapping 
works by sending the point at z = to u — > — oo, thereby 
folding the line at Im z = in half so as to form the two 
edges of the strip. Thus, the boundary condition of no 
supercurrent through the upper and lower edges of the 



strip, i.e., dh a {u)4'{ u )\ u _{ w y becomes the condition of 
no current through the line Imz = 0, i.e., d\ m (z)4>{ z ) | _ - 
The phase field (f>(z; zq) = Im [ln(z — zq) — ln(z — zjj)] of 
a vortex located at z = xq + iy§ in the upper half-plane 
can be found using a single image- vortex. Next, by us- 
ing the Cauchy-Riemann equations, we identify the vor- 
tex potential Ghaif-pianc as the real part of the harmonic 
function whose imaginary part corresponds to the vortex 
phase field, i.e., 

Ghaif- P ianc(z; z ) = 2np s Re [ln(z - z )-ln(z - Zq)] . (3) 

Applying the conformal mapping, we find the result 



Gstripfu; w ) = 

27tp s Re [\n{e™ /w 



(4) 



a TTU /w\ 



ln(e 



TTu/w TTUq/w\ 



for the potential G str i P of a vortex in the strip. Using this 
vortex potential, we find that the energy due to inter- 
vortex interactions is given by 



E r 



2 



Gstrip( u i', u j)- 



(5) 



We note that the inter-vortex interaction becomes zero 
for vortices close to the top or bottom edges of the strip. 
Thus, the vortices in the strip do not influence the po- 
tentials of the vortex reservoirs. 

Vortex blockade - We begin by examining the Weber 
(cf. Coulomb) blockade in the absence of a bias current 
(cf. source-drain voltage). The key to either blockade is 
that for generic values of the magnetic field B (cf. gate 
voltage V g ) the energies for N or N + 1 vortices (cf. elec- 
trons) to be on the strip (dot) differ, and therefore vortex 
(cf. charge) current does not flow. However, for special 
values of B (cf. V„), the degeneracy condition 



E N (B,J = Q)=E N+1 (B,J = Q) 



(6) 



is met, allowing vortices (cf. electrons) to move freely on 
to and off of the strip (cf. dot), thus lifting the blockade. 

For the case of nonzero bias current J, we must con- 
sider separately the processes of adding vortices from the 
top and bottom reservoirs. These considerations furnish 
two stability conditions, viz., 



En(B, J) + Stop 

En{B, J) + -^bottom 



En+i(B, J), 
En+i(B, J), 



(7a) 
(7b) 



that define the Weber diamonds, in which E top and 
-^bottom are the potentials of the top and bottom reser- 
voirs. In Fig. [3] we plot the first few diamonds in the 
J vs. B plane. The regions inside the diamonds corre- 
spond to 1,2,... vortices on the strip in the blockaded 
regime (i.e., poor vortex conduction across, and there- 
fore good electrical conduction along, the strip). For in- 
creased bias currents, the potential difference between 
the vortex reservoirs overcomes the blockade, driving the 
strip into the good vortex conduction regime (i.e., to poor 
electrical conduction along the strip). 
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FIG. 3: Weber diamonds in the bias current vs. magnetic 
field plane (w = 94 nm, L — 1 /^m, £ = 18 nm, and p = 5K). 
Regions inside the diamonds correspond to small electrical 
resistance along the superconducting strip or, equivalently, to 
small vortex conductance across it). 



reservoirs but in the strip itself. However, there are many 
other possible damping mechanisms, such as coupling to 
the bath of Bogoliubov quasiparticles. The question of 
which mechanism dominates depends on the microscopic 
properties of the superconductor. We shall focus here on 
the simplest possibility, which is to use as the thermal 
bath the normal modes of the line of vortices inside the 
strip. 

In our phenomenological approach, we begin by mod- 
eling the line of vortices inside the strip as a Luttinger 
liquid having a single gapless mode, corresponding to the 
longitudinal sound mode of the one-dimensional vortex 
crystal. Next, we describe the coupling between the vor- 
tex reservoirs and the Luttinger liquid by introducing a 
model Hamiltonian that accounts for vortex tunneling 
between them: 



In order to estimate the actual magnitude of the con- 
ductance of the superconducting strip, we need to extend 
the static model of the blockade to include the dynam- 
ics of vortex hopping and energy dissipation. For the 
quantum-dot case, dissipation arises in the leads, which 
act as "decohering baths" in which the quasiparticles are 
assumed always to obey the Fermi-Dirac distribution [5] . 
In the vortex case, dissipation occurs not in the vortex 



H tun = J dx b(x) (t t e l6t + t b e t9b ) + h.c. , (8) 

where b(x) is the operator that annihilates a vortex 
at position x along the strip. Following Beenakker's 
approach to quantum dots [5J, which invokes Fermi's 
Golden Rule, we obtain a master equation for the proba- 
bilities Pq, Pi, ■ ■ ■ of having 0, 1, . . . vortices on the strip: 



d t P N = PN+ a e~ Pei U(e/|#tunhW| 2 T a S[(E N+a + ei)-(E N +ef)-a£ a ] (9a) 

e il e/cr=±l aG{t,b} 

- J2 P NC~' 3 ^ |jv+ ff (e/|tftun|e*)w| 2 ^ a S[(E N + ei)-(E N+a + e f )-a£ a ] (9b) 

Ci,€f cr=±l ae{t,b} 



where \ei) N and \ef) N correspond to the states of the 
Luttinger liquid (with N vortices), £ t and £ b are the po- 
tentials of the top and bottom vortex reservoirs, and the 
(5-functions account for energy conservation. 

In the regime of small bias-supercurrent J (cf. low 
source-drain bias voltage) , we can solve the master equa- 
tion and use the solution to recover the vortex conduc- 
tivity a (cf. electrical resistivity): 



a{J) 



Tt+T b 



dte lJt 



dx dx' 



(t 2 + (x- X 1 ) 2 ) 



(10) 



where K is the corresponding Luttinger parameter. Here, 
we have used the standard trick |13) of converting the 
(5-functions into time integrals to recover the retarded 
Green functions. In general, the [x x)' integrals will be 
cut off by the system size at large argument and by the 
inter- vortex spacing at small argument. Finally, we com- 
ment that the boundary conditions for the Luttinger liq- 
uid are determined by the coupling of the vortices in the 
strip to the degrees of freedom in the leads. 



Comparison with experiment - Using Eq. ([6| , we can es- 
timate theoretically the period of the magnetoresistance 
oscillations for the sample geometries of Ref. [7j. The 
only undetermined parameter in the theoretical estimate 
is the coherence length. Treating the coherence length as 
a fitting parameter, we were able to fit the periods of all 
four narrow-strip samples (Table [l]) . We remark that the 
experimentally observed period is always longer than the 
geometric period associated with adding a flux quantum 
to the area of the strip (Table [l|. Our theory implies that 
the effective width of the strip is narrower than its geo- 
metrical width, due to the strong attraction of vortices 
to the strip edges by their images (i.e., Bi mage ), thus re- 
solving the period puzzle. The experimental signal shows 
only mild magnetoresistance oscillations, rather than the 
pronounced vortex-blockade peaks that we predict. We 
suspect that this smearing is a result of T t and r fc being 
large, and thus the strip operating in between the open- 
and closed-dot regimes. Alternatively, the smearing may 
be a result of thermal broadening. If so, by lowering 
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TABLE I: Coherence length £ (a fitting parameter) required 
to reproduce the experimentally measured period P using 
our theoretical model for four samples [12] of widths w and 
lengths L (Table I of Ref. [7]). In each case, we fit the ex- 
perimental period at the intermediate value of magnetic field 
B ~ 1.5 T. For comparison, we also provide the geometric 
period P gC0 mctric = (&owL)~ . For each sample, we find a 
reasonable value of coherence length, which should lie in the 
range from 10 nm to 30 nm, according to Ref. [7] . 



the temperature the blockade features are expected to 
be strengthened. 

Concluding remarks - We have demonstrated that mag- 
netoresistance oscillations in superconducting strips can 
be readily be described using vortex — as opposed to 
Cooper-pair or charge — coordinates. Via an analogy 
with the physics of quantum dots, we have constructed a 



model for the 'Weber blockade' of a superconducting strip 
that captures the essential features of the magnetore- 
sistance oscillations observed experimentally. The exact 
dissipation mechanism remains to be understood; it will 
determine the shape of the Coloumb blockade peaks pre- 
dicted in Eq. ( 10 ) but not their locations. Alternatively, 



one could use the experimental blocakde signatures in 
order to investigate the nature of dissipation, according 
to Eq. ((TO I. 
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